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Abstract 

The correlation between the values of wavefunctions at two different spatial 
points is examined for chaotic systems with time-reversal symmetry. Em- 
ploying a supermatrix method, we find that there exist long-range Friedel 
oscillations of the wave function density for a given eigenstate, although the 
background wavefunction density fluctuates strongly. We show that for large 
fluctuations, once the value of the wave function at one point is known, its 
spatial dependence becomes highly predictable for increasingly large space 
around this point. These results are compared with the experimental wave 
functions obtained from billiard-shaped microwave cavities and very good 

agreement is demonstrated. 
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Quantum properties of classically chaotic systems such as "billiards" and quantum dots 
are revealed to have remarkable universal behaviors which depend only on the generic sym- 
metry of the system, such as the time-reversal symmetry and/or the spin rotational sym- 
metry It has been shown that the spectral statistics are well described by universal 
statistical correlations derived from the random matrix theory (See also for re- 
views.) Complementary and comprehensive information beyond the energy statistics can 
be obtained by examining the statistics of chaotic wavefunctions. For example, the distri- 
bution of the local density in a fully chaotic system is known to be universal and obey the 
Porter-Thomas distribution which is given for a system with time-reversal symmetry, 
by the equation 

Poiv) ^ {5iv - V\M^)n = ^ exp(-^;/2), (1) 

where V'el'') is the eigenfunction with energy e in a system with volume V. (■ ■ ■) means 
the average over the disorder and/or irregular potential. Eq. (0) tells us that wavefunctions 
fluctuate strongly but in a universal way. 

To get further understanding about the nature of chaotic wavefunctions, other statistical 
quantities which can characterize their spatial correlations are desirable and needed. The 
average behavior of the amplitude of the wavefunction has been conjectured by Berry to be 
similar to a pattern generated from random super-position of plane waves [|lO],|ri|. Based 
on this assumption the average amplitude correlations were shown to be a Bessel function. 



Recently, these expressions have been derived within the supersymmetry formalism |]12 
In this paper, taking a chaotic system with the time-reversal symmetry, correlations about 
a particular value of the wavefunction have been derived analytically and compared to 
experiments for the first time. This is not only a more stringent test of the universality of 
chaotic wavefunctions, but also gives us a handle on knowing the behavior of a wavefunction 
of the system, once the wavefunction is known only at a limited number of points. Our 
main object is to investigate the joint probability distribution function of the density for 
two different spatial points (r = |ri — r2|) defined by 



(2) 

Although the relevant universality class for experiments which can directly observe the 
amplitude of wavefunctions, such as quantum corrals ITslJT^ and microwave cavity HTsjpiB 



IS 



orthogonal, only the expression for P{vi, V2) in the unitary case is known so far ||T^, because 
of technical difficulties. Here we evaluate P{vi,V2;r) for the orthogonal case by finding 
special techniques (Eq. (|20|) below). In the microwave cavity, the electromagnetic field 
obeys the same equation of motion as a quantum particle in a two-dimensional billiard. This 
enables us to make direct comparison between the analytical results and the experimental 
data. The experimental data for the wave function density was obtained from thin cylindrical 
microwave cavities of the Sinai stadium, by using a cavity perturbation technique |16|,|T8[. 
The wavefunction density data was earlier seen to be consistent with Eq. (|l]) and also 



in agreement with the expression for density auto-correlations obtained in ref. [0. Once 
we know P{vi,V2',r), we can also find the conditional probability 

P,,{v2;r)=P{vuV2;r)/Po{v,). (3) 

Pvi{v2;r) describes the distribution of the wave function of V2 = y|V^e(i"2)P at the point r2, 
provided that Vi = V\ipe{'ri)\'^ at ri. Here we compare the coordinate dependence of the first 
and the second moments of Eq. @ between theory and experiment. 

The system under consideration can be expressed by the Hamiltonian 

i7=-^p2 + f/o(r) + f/i(r), (4) 
2m 

where Uo{r) denotes the regular part of a confining potential, and Ui{r) is a random po- 
tential which is responsible for the chaotic dynamics — impurities or 'imperfection' of the 
shape of the system. We take the ensemble average over f/i(r) by use of the supermatrix 



method, which reproduces the spectral correlations of Wigner-Dyson statistics p!9| , |20| and 
recently was successfully applied to calculate other universal properties relating to chaotic 
wavefunctions OJT7|J2T|J2! 



We should remark that although the supermatrix method was originally derived from 
the Gaussian random potential where the mean free path £ is much smaller than the size of 
the system L, the ergodicity hypothesis allows us to extend our present result for cases 
£ ~ L by identifying the averaging over space and different states for a given sample with 
that over disorder. As a confirmation of the ergodicity hypothesis we will demonstrate that 
the theoretical dependencies for wavefunctions derived from a disordered system with £ <^ L 
are universal and describe very well experimental results for quantum billiard systems for 
which £ ~ L. 

After lengthy calculations which we will sketch later, we have obtained the following 
analytical expression for P{vi,V2',r) for the orthogonal case: 



1 f/M pdp d\ l-p^ p dz ,2 ^ ^i + z V2 + Z 

2nf{r) Jo f2(r\-p2 \ dp) \ 2'kviV2 h 2 2 



where /o(p) is the modified Bessel function and /(r) is the Friedel function ||2^. Note 



that P(fi,f2;r) depends on /(r) in an universal way and all coordinate dependence of 
P(i;i,f2;r) is incorporated only through fir). In fact, the function fir) represents the 
average correlation of the amplitude of wave function |TO| , p!Tl| : /(r) = V^ < il'ei''^i)'^e{f^2) >■ 



For the case of a flat background potential in a (i-dimensional system, /(r) becomes 

/(r) = T{d/2) J„2-i{kr) e-^'^\ (7) 

where k is the wave vector (e = h^k"^ /2m), Jn{x) is the Bessel function, and r(n) is the 
gamma function. Note that the envelope of /(r) decays like (fcr)^^'^"^^/^ for k~^ ^ r ^ i 
and this behavior corresponds to the representation of chaotic wave function as a random 
superposition of plain waves pi]|JTT| . 

It should be remarked that the same distribution function in the unitary case is given 



by P{vi,V2;r) = (p{vi,V2; f{r)) [T^. Due to two additional integrations in Eq. (H) the 



4 



spatial correlations for the orthogonal case are weaker and the fluctuations are stronger in 
comparison with the unitary one. 

We flrst check that P{vi,V2',r) given by Eq. yields correct limiting behaviors. For 
remotely separate points such that /(r) ^ the fluctuations of the wavefunction density 
within the given eigenstate become independent, i.e., P{vi,V2;r) ^ Po{vi)Po{v2) ■ In the 
opposite limit of close enough points that /(r) ^ 1, there is an obvious strong correlation 
between fluctuations as 



P[Vi,V2;r) ^ , = exp 



(8) 



From Eq. (|^), we can extract information about the gradient of the wavefunction. By 
setting 4'{r2) = '?/'(ri) +rVn'?A(ri) + C'(r^), and expanding for small r, we obtain the joint dis- 
tribution involving the wavefunction and its gradient in any direction n = r/r. Accordingly 
we flnd that the gradient of the wavefunction along any direction fluctuates independently 
of the value of wavefunction, and obeys also the Porter-Thomas distribution: 

{6{v - Vmr)ms - ^iVnV^WH) = Po{v)Po{s). (9) 

This conclusion is, however, not true for higher gradients of the wavefunction. 

The conditional probability P^^(f2;r) is obtained straightforwardly with Eqs. To 
see how the fluctuations of the wavefunctions behave and to compare between the analytical 
results and the experiments, the conditional average (^2)1,1 and the conditional variance is 
more convenient. Denoting 6v2 = f 2 — (^2)1,1, we obtain: 

{v2)., = l + f{r)iv^-l), (10) 
{{Sv2f)., = 2 + 4f (r) {v^ - 1) + 2f\r) (1 - 2v,). (11) 



Comparing with results obtained for unitary case we flnd that the conditional average 
Eq. (p!0| ) is exactly the same. Thus we cannot tell the symmetry of the system only from the 
averaged amplitude even if we know the conditional one. To detect the sjTumetry, we have 
to examine the variance, where there is a factor of 2 difference between the orthogonal and 
the unitary case. 



In Fig. and 0, we compare the analytical results of the conditional average and variance 
with the experimental data from microwave cavities. The experimental curves were obtained 
by picking points in a wavefunction with the same value and calculating the average wave- 
function value a distance r from it on a circle. This quantity was then again averaged over 
at least 50 wavefunctions after rescaling the wave number to obtain better statistics. Fig. |I| 
shows the plot of {v2)v^ in Eq- (0) with experimental results for vi = 2 and vi = 7. Very 
good agreement is seen for both sets. In fact, agreement is excellent for all values of vi above 
1, below which the noise and errors in the measurement of the wavefunction measurements 
lead to qualitative differences. In Fig. |^, the comparison of the data to the expression in 
Eq. (|lT]) was done. Again one sees an excellent agreement with experimental errors of 5 %, 
which is the level of experimental accuracy. 

According to Eqs. (|10|,|lT]), we can say, as in the unitary case that large fluctuations 



of the wavefunction have some striking structure which is not present for small fluctuations. 
For f 1 ^ 1 the ratio of the variance to the average square is 

^ - - n . (12) 

Therefore at r ^ ^, where the "correlation length" C, ~ k^^vl^^'^ ^ k~^, the variance 
{{Sv2)^)vi can be very small in comparison with (^2)^^- It means, that once we know that 
the wavefunction |^/'(r)p is equal to Vi at ri, it is highly likely to have a value (f 2)1^1 ~ pvi 
at r2 for r ^ ^. In this sense, the large fluctuation behavior of the wavefunction becomes 
highly predictable. In contrast, for small fluctuations V2 <^ 1, we easily see that {{Sv2y)vi ~ 
2(^2)^^, independent of vi. We also find directly from Eq. (|]) that fluctuations turn out 
to be independent, i.e., P„^(f2;r) ^ Although more careful evaluation gives us a 

correlation length ^ which ensures independent fluctuations for the region r ^ ^, ,^ in such 
an evaluation turns out to be very small, i.e., C, ~ k^^^/ih ^ k"^. 

These behaviors are qualitatively the same both in systems with or without time-reversal 
symmetry. Therefore we can say that this is a generic property of chaotic wavefunctions. 
Also, in the semiclassical description of chaotic systems, periodic and closed orbits are known 
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to be associated with large values of the wavefunctions |'?/'(r)p [^,^. In this respect, our 



present results may imply there is some structure present in these orbits. 

Now let us proceed with the derivation of our main result in Eqs. To evaluate the 

joint distribution P{vi, V2',r), we work with its moments, 

qnmir) = V^^^ (|V'.(ri)p"|V'.(r2)|''") = (^^W)- (13) 

qnm is known to be closely related to the moments of the exact retarded and advanced Green 
functions G^'^, 



(14) 



where v is the average DOS and 



We can obtain qnm in terms of Fnm by the relation [25,26 



- 1)! (m - 1)! , /7\"+'"-i 

i^. X Fnm, (16) 



2(n + m-2)! 7-0 V A. 
since the leading contribution to Fnm. for small 7 comes from the state whose energy ea 
coincides with e. 

Fnm can be evaluated by the supermatrix method. However, since we cannot utilize a 
simple expression like Eq. (12) of Ref. [jl^ for the orthogonal case, we are forced to expand 



Fnm directly by the Friedel function /(r) as (see also Ref. |lT2tpT 



Fnm{r]l)=n\m\Y, C,{n,m) f^ir). (17) 



Using the same notation for the supermatrix elements as in Ref. [|T^, and defining Q = 
(^_lYQab (fQ^ a = 1, 2), the coefficient C,( n, m) in Eq. ( p^Tf ) is given in terms of Q-supermatrix 
elements: 

C„(n,m)=V n 7^ F^r-r TH VT^ (18) 



where the summation is taken over the all possible combinations of nonnegative integer kab, 
Pab, Lb which satisfy the condition: 2q = Y.a^b{Kb+Pab + lab), m-Y.aha = '^hi + Y.a^bkab = 
2pii + J2a^bPab, and n - Ea^a2 = 2A;22 + Y^a^bkab = 2p22 + Ha^bPab- The symbol (■ ■ ■)q 
denotes an integration over the saddle point manifold, i.e., 

(■ ■ .)q ^ JdQ (■ ■ •) exp[-(7r7/4A)StrAg], (19) 



where the definitions of A and Str as well as the structure of the Q matrix are found in |]T9 

In principle, averaging (■ ■ ■)q in Eq. (0) can be carried out by using the parameterization 
of Ref. [|19|. However, we found it technically unfeasible to evaluate this expression in such a 



general form. Fortunately, to get Eq. (|1^) one need only to know Fnm{l 0). The leading 



contribution in this limit can be extracted by transforming parameters Aj = l + MiyA/7 for 
2 = 1, 2 in the parameterization given in Ref. [1^. By calculating the leading order of small 



7, we find the relation 

Qcd Qc'd' — Qcd' Qcd'J (20) 

for a,b = 1,2 and c,d = 3,4. Substituting Eq. ( pUD into Eq. (|T^ and combining with 
Eq. ([TBp, we finally obtain after integrating over the Q matrix, 

{2n-l)\\{2m-l)\\p'>{r) 

Qnm[r) 2^ ^n+m-2q _ ^)! _ g), (2g)! ' ^^^^ 

Reconstructing P{vi,V2',r) from the moments qnm completes the derivation of our main 
result Eqs. (|,|). 

In conclusion we have presented analytical results for universal statistical quantities 
which characterize the coordinate dependence of chaotic wavefunctions of the system with 
time-reversal symmetry. Further we have demonstrated excellent agreement between the 
theoretical results and experimental results of microwave cavities. The spatial correlations 
demonstrate the long-range Friedel oscillations of wavefunction density and the existence of 
extended spatial regions of high wavefunction density. 
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FIGURES 

FIG. 1. The average spatial dependence of chaotic wavefunction squared {v2)vi = I^(|^e(r2)P) 
whose value at a point ri is known to be f i ( r = |r2 — ri|, V \s the volume, and 

h'^k'^ / {2m) = e). The theoretical prediction Eq. ( p!o| ) is compared with experiments from the 
microwave Sinai stadium cavity for vi = ?(□) and vi = 2> (o). Inset: Representative eigenfunction 
of the chaotic Sinai stadium. 

FIG. 2. Comparison of the conditional variance of a wave function {{5v2f')v-i {^^2 = V2 — {v2)vi 
and V2 = V^|V'e(r2)P) as a function of the distance from point ri and a reference value vi 
{vi = y|'(/'e(ri)P) between the theory (Eq. (|Tl])) and the experiment for vi = 1 (□) and f i = 3 (o). 
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